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Mixed order parameters, accidental nodes and broken time reversal symmetry in 
organic superconductors: a group theoretical analysis 

B. J. PowelQ 

Department of Physics, University of Queensland, Brisbane, Queensland 4 072, Australia 

We present a group theoretic analysis of several classes of organic superconductors. We argue 
that highly frustrated half-filled layered organic superconductors, such as k-(ET)2Cu2(CN)3 (where 
ET is BEDT-TTF) and /3'-[Pd(dmit)2]2-'f , undergo two superconducting phase transitions, the first 
from the normal state to a d-wave superconductor and the second to a d + id state. We show that 
the monoclinic distortion of k-(ET)2Cu(NCS)2 means that the symmetry of its superconducting 
order parameter is different from that of orthorhombic «:-(ET)2Cu[N(CN)2]Br. We propose that 
quarter filled layered organic superconductors, e.g., S-(ET)2l3, have dxy + s order parameters. 

PACS numbers; 74.20.Rp, 74.70.Kn, 02.20.-a 
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One of the most basic questions that can be asked 
about any phase of matter is, what symmetries does 
it spontaneously break? For example, all superconduc- 
tors break gauge symmetry but many also break addi- 
tional symmetries 0, Q- Organic charge transfer salts 
are an important class of superconductor because they 
are highly tunable and have a number of exotic proper- 
ties WM ^, "6] such as a small superfluid stiffness |y,l3j 
a Mott transition spin liquid states l E , IQ, llj , a 'bad 
metal' , pseudogap like behaviours |l2l |. and charge or- 
dered states 01 . Many of these properties are analogous 
to the cuprates Q| and this has heightened interest in 
organic superconductors as model systems. Yet, in spite 
of the two decades of intense effort to understand these 
systems, many basic questions remain unanswered. Of 
particular importance is the, still controversial, question 
of what is the symmetry of the superconducting state 

Ellii. 

In this Letter we use the methods of group theory to 
analyse the symmetries of the superconducting states in 
these materials. We discuss the pairing symmetries of k- 
(ET)2Cu[N(CN)2]Br (henceforth K-Br) and k-(ET)2Cu- 
[N(CN)2]C1 (k-C1) which we argue are 'd.^2_y2\ and k- 
(ET)2Cu(NCS)2 (k-NCS) which we show has additional 
components of other 'd-wave' states and thus accidental 
nodes (i.e., nodes that are not required by symmetry). 
Extending these arguments allows us to show that highly 
frustrated materials such as k-(ET)2Cu2(CN)3 (K-CN3) 
and /3'-[Pd(dmit)2]2-'^ (/3'-dmit) [ij wiU undergo two su- 
perconducting transitions the first from a normal metal 
to a 'd-wave' and the second from a 'd-wave' supercon- 
ductor to a 'd + id' state which breaks time reversal sym- 
metry (TRS). Similar reasoning implies that the quarter 
filled organic superconductors have ^dxy + s' order pa- 
rameters. 

Group theory provides a powerful tool for addressing 
the symmetries of superconducting states as it does not 
assume any particular microscopic mechanism or theory 
of superconductivity P, Q| ■ Such approaches are vital for 
the organic superconductors where there are clear signs 



that BCS theory, and indeed weak coupling approaches in 
general, are not sufficient to explain the observed experi- 
mental results, particularly the small superfluid stiffness 
0, 0- Further, the properties of nodal quasiparticles, 
which are determined by the symmetry of the supercon- 
ducting state, have proved crucial in both the supercon- 
ducting and pseudogap states of the cuprates |l8| . 

Models of superconductivity based on both electron- 
phonon coupling |l5| and the Hubbard model on an 
anisotropic triangular lattice 0, have been pro- 
posed for the half-filled layered organic superconductors 
(^LOS). Calculations based on these two microscopic 
models lead to different conclusions. For phononic mech- 
anisms both isotropic, nodeless 's-wave' states belonging 
to the Ai irreducible representation (irrep) of C2V and 
'da;2_y2' states belonging to B2 with nodes in the gap 
along the lines kx — ky and kx ~ —ky have been pre- 
dicted. Calculations based on the Hubbard model sug- 
gest a (B2) 'd^2_j^2' state |5lllll9t. 

Experimentally, the pairing symmetries of K-Br, k-C1, 
and K-NCS are unclear. All show signs of unconven- 
tional superconductivity: there are no Hebel-Schlicter 
peaks jl2j , the thermodynamic measurements performed 
to the lowest temperatures show a power law tempera- 
ture dependence 20] and the disorder strongl y su ppresses 
the superconducting critical temperature |l6ll21| . Triplet 
superconductivity can be ruled out on the basis of 
measurements of the Knight shift |l2 1 and upper critical 
field. 
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^LOS. Both K-Br and k-C1 are 
with orthorhombic (-D2/1) crystal structures. There are 
four even parity irreps of D2h (see table IJ) all of which 
are one dimensional. Canonically, the highly conducting 
plane is the ac plane in both K-Br and k-C\. Further, 
it is usual to define the x- and y-axes as lying along the 
directions of the largest inter-dimer hopping integrals, 
which lie along the diagonals of the ab plane so that x = 
(b + c)/2 andy = (b - c)/2. 

^LOS crystals are extremely anisotropic: the inter- 
plane hoping integral is about three orders of magni- 



Irrep Required nodes Example basis functions States Example basis functions (fcz||c) States 

Aig none Ik, Ai, B'i, C^, XkFk s, d^y Ik, Xkik, Ik + -fkik s, d^^y 

Big line j4k-Bk, (^k + ik)^k d(x+y)z ^kSk, Xk — Fk d^^-y^ 

B2g line AkCk, Xk — Fk d^2_y2 ^kCk, (Xk + lk)^k d(a;+s;)z 

B3g line B^C^, (Xk — Y]^)Zi^ rf(3:-y)z j?kC'k, (Xk — Yi^)Zi^ d(^^_y)^ 

TABLE I; The symmetry required nodes of the even parity irreps of the point group D2h, which represents the symmetry of 
the orthorhombic organic superconductors such as K-Br, in which the highly conducting plane is the a-c plane, i.e., kz\\b, and 
S-I3, in which the highly conducting plane is the a-b plane, i.e., k^\\c. The functions Ik, Xk, Kk, -^k, ^k, Bk and Ck may 
be any functions which transform, respectively, as 1, k^, ky, kz, ka, kt and kc under the operations of the group and satisfy 
translational symmetry. 



tude smaller than that in-plane [Sj . The superconducting 
properties are also extremely anisotropic Hence, su- 
perconductivity in which the order parameter transforms 
as either the Big or B^g irreps are extremely unlikely Q . 
Therefore, a ^(1^2 _y2^ state transforming as the B2g ir- 
rep of £>2/i is most consistent with the experimental and 
theoretical evidence. 

Monoclinic ■^LOS. The best studied ^LOS with a mon- 
oclinic crystal structure is k-NCS. This material has a 
C2/1 point group and the highly conducting plane is the 
be plane. Here, the x- and y-axes are usually taken to 
be along the diagonals of the 6c-plane. The only non- 
identity irrep of C2/1 that corresponds to singlet super- 
conductivity is Bg (see table . This has symmetry 
required nodes only along the c-axis. 

However, a subset of the possible choices for the basis 
function of the Bg irrep lead to ^d^2_y2^ superconductiv- 
ity (i.e., states with nodes along fc^ = fc^), which we have 
just argued is the superconducting state realised in the 
orthorhombic ^LOS. Further, the properties of k-NCS 
are so similar to those of K-Br that the differences be- 
tween the two materials are often described as 'chemical 
pressure' [El. Therefore, one expects that the supercon- 
ducting states of the two materials are closely related. If 
the node is shifted away from the 6-axis in the higher sym- 
metry, orthorhombic, case (as sketched in Fig. ^ then 
there is a change in symmetry which is accompanied by 
a phase transition. There is no requirement for such as 
phase transition in the monoclinic case as both possible 
the order parameters shown in Fig. ^ transform as the 
Bg irrep of C2h ■ Any finite contribution to the supercon- 
ducting order parameter from basis functions which do 
not have nodes along 6-axis will cause this node to move 
or, for a sufficiently large contribution, disappear. Thus, 
predictions or measurements which show that monoclinic 
^LOS have nodes anywhere except along the c-axis are 
not robust. Therefore one cannot reliably conclude that 
the 'd-wave' {Bg) state is ^d,j.2_y2\ If there is an acci- 
dental node then it will not lie along 6-axis and it may 
have either a temperature dependence, a dependence 
or both. 

The above argument also shows that calculations based 
on the anisotropic triangular lattice (whose symmetry is 
represented by the group C21,) excludes components of 
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FIG. 1: (Color online.) Sketch of the differences between the 
nodal structures of (a) orthorhombic ^LOS (e.g., K-Br) and 
(b) monoclinic ^LOS (e.g., k-NCS). In both cases the Fermi 
surfaces are shown as solid lines and the nodes in the gap 
are shown as dashed lines. For d^2_y2 states is orthorhombic 
^LOS symmetry requires that the nodes lie along the a- and 
c-axes. Whereas for d-wave states in monoclinic ^LOS the 
only required symmetry is that the order parameter be anti- 
symmetric under rotation by tt about the 2-axis. As a node is 
not required along the 6-axis in the 'd-wave' state of the mon- 
oclinic iLOS we expect contributions from basis functions 
that do not go to zero along the fo-axis to move this node and 
possibly remove this node altogether. High symmetry points 
are given their canonical labels. 



the gap in a manner that is not relevant to the monoclinic 
^LOS. Thus weak interactions excluded by these models 
will dramatically change the symmetry of the gap. Even 
neglecting such interactions, at sufficiently low temper- 
atures at least small additional 'c?(a;-j,)z' components to 
the order parameter should be expected. Exotic states 
with more nodes than are present in '(i^2_j^2' states have 
been proposed 22]. Such nodes would almost certainly 
be lifted in both monoclinic and orthorhombic ^LOS as 
they are far from robust. 

TRS breaking in K-CN3 and (3' -dmit. Much atten- 
tion has been focused on K-CN3 following the dis- 
covery that, in spite of there existing well formed local 
moments, in the low pressure, insulating phase, these 
moments do not order down to the lowest temperatures 
studied (20 mK) 0. Both Huckel calculations "231 and 
fits of the susceptibility calculated by series expansions 
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line 



^kCk, (Xk - Yk)^k, 
-BkCk, — Yk 



TABLE II: The symmetry required nodes of the even-parity 
irreps of the group C2h which represents the symmetry of n- 
NCS and several other charge transfer salts with monoclinic 
unit cells. Note that the symmetry line node in the Bg irrep 



is required to lie in the plane fee 
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0, 13 to that observed experimentally suggest that the 
band structure of K-CN3 is that of the isotropic triangu- 
lar lattice. Series expansions 0, 0| also show that the 
triangular lattice Heisenberg model is a good approxima- 
tion for the low pressure, insulating phase of /?'-dmit. 

The symmetry of the triangular, or more correctly 
hexagonal, lattice is represented by the C^v point group. 
An interesting feature of C^y is that it has two two- 
dimensional {2D) irreps. In a 2D irrep the order param- 
eter, Ak, is a linear combination of the basis functions, 
^-k' , of the irrep, i.e., Ak = Vi'^k + ^2^k- Hence, on 
the hexagonal lattice the Ginzburg-Landau free ene_rgy 
of order parameters belonging to the 2D irreps is 



energj 

m 



+(32{v*im - mv2f- (1) 

The ground state solution, ff = (771,772), is: (i) ?7 = (1,0) 
or (ii) ff = (0, 1) for P2 > (the degeneracy is lifted by 
sixth order terms 1, 2]); (iii) 77 — (1,?) for /?2 < (this is 
the weak coupling solution). In the E2 irrep ^k describes 
a 'dj.2_y2' state and ^k describes a 'dxy^ state. Thus the 
three solutions correspond to (i) 'dx2_y2' superconductiv- 
ity; (ii) 'c?xi/' superconductivity; and (iii) 'd^^-y-^ + idxy 
superconductivity. Thus state (i) is the same state as we 
have discussed above for the orthorhombic ^LOS. How- 
ever, a number of studies [2^ of the Hubbard model on 
the hexagonal lattice suggest that state (iii) is realised. 

The role of monoclinicity. It is important to realise 
that although the band structures of K-CN3 and /3'-dmit 
are close to the hexagonal lattice, they will not be pre- 
cisely that of the hexagonal lattice because these materi- 
als form monoclinic crystals. The monoclinic distortion 
of the crystal lowers the symmetry of the microscopic 
Hamiltonian. We account for this perturbation by in- 
troducing a symmetry breaking field, e which, to lowest 



order, enters the free energy as £(|77i| 



\V2\ 



2) "2^. 



Thus, 



a+|77i|2 + a_|r72|2+/3i(|77i|2 
+P2{rilV2 -ViV2f, 



2\2 



\mV) 



(2) 



where a± = a{T — Tc±) and Tc± — Tci: e/a. This the- 
ory predicts that the monoclinic crystal will have two 
superconducting transitions, the first to either a 'dx2_y2^ 
superconducting state or a ^dxy' superconducting state 
(both of which have nodes) and the second to a, fully 
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FIG. 2: Sketch of the proposed phase diagram for supercon- 
ductivity near the hexagonal lattice in K-CN3 and /3'-dmit. e 
is a symmetry breaking parameter which lowers the symmetry 
from C^v to C2h- Physically e could represent uniaxial strain 
or pressure, e 7^ at ambient pressure due to the monoclinic 
crystal structure. 



gapped, 'dx2-y2 -\- id^y state. This leads us to propose 
the phase diagram sketched in Fig. |21 For small e the 
difference in the two Tc's grows linearly with e. But, for 
large £ the symmetry of the lattice will cease to be related 
to that of C^y and return to that of C2V Therefore, for 
large e all of the irreps are one-dimensional and, neglect- 
ing the possibility of accidental degeneracies, we expect 
a single superconducting transition. A similar scenario 
has previously been studied in some detail in the con- 
text of the double superconducting transition observed 
in UPta l^^. However, for UPta the proposed symme- 
try breaking field arises from a weak antiferromagnetic 
background 0] , rather than from the crystal structure. 

The above predictions are readily testable. Any num- 
ber of experiments might see the double superconducting 
transition (e.g., specific heat or ultrasound). Further, the 
proposed low temperature 'd + id' state breaks TRS, this 
could be detected directly by a number of experiments 
1] most notably /iSR |27l |. Experimental confirmation 
of a double superconducting transition and broken TRS 
would be extremely important, not only because of the 
intrinsic interest in these phenomena, but also because 
they would be conclusive proof of unconventional su- 
perconductivity in the layered organic superconductors. 
The implication is then that phonons are not responsible 
for the superconductivity in these materials, but rather 
strong electronic correlations. 

Quarter filled materials have been less studied ex- 
perimentally, but, measurements of the in- and inter- 
plane penetration depths show power law dependencies 
suggesting unconventional superconductivity |28l| . Some 
quarter filled layered organic superconductors (|LOS), 
e.g., 0-(ET)2l3 {O-I3), have unit cells with a D2h point 
group, while others, e.g., /3"-(ET)2SF5CH2CF2S03 (^- 
SO3), have d point groups. It has been proposed [iJl 
that a minimal model for |:LOS is the quarter filled ex- 
tended Hubbard (or t-U-V-J') model on the square lat- 
tice. (Note that the unit cell used in this model is ro- 
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Point group Example material 


Irrep 


State 


C2h {Cqv) 




B, + iA, {E2) 


'd + id' 


D2h 


K-Br 


-B29 




C2h 


K-NCS 




^dj.2_y2 + (i' 




e-h 




dxy ~\~ S 


a 


/3"-S03 


A, 





TABLE III: Summary of the superconducting states proposed 
on the basis of the group theoretic analysis in this Letter. The 
parenthetic point group and irrep in the first row indicate the 
approximate symmetries which drive the physics. 



tated with respect to the crystallographic unit cell so that 
X — (a+b)/2 and y = (a— b)/2.) This has led to the idea 
that charge and spin fluctuations may cooperatively me- 
diate ^dxy superconductivity However, this predic- 
tion of the location of the nodes is not robust as the crys- 
tal lattices have significantly lower symmetries than the 
model. In the D2h point group (tabled} the basis func- 
tions that describe the 'd^i/' superconducting state belong 
to the trivial Aig irrep. Therefore, these nodes are not 
robust in the same sense as 'd^2_j,2' states are not robust 
in monoclinic ^LOS . Thus one expects that the ^LOS 
will have an ^d^y + s' order parameter. The order pa- 
rameter may still have accidental nodes, but such nodes 
will be shifted away from the lines k^ky = 0. The same 
arguments and conclusions hold for -jLOS with Ci point 
groups. A similar argument suggested that cuprates with 
a small orthorhombic distortion are 's -f dxi_y2^ super- 
conductors '2^ . Tunnelling experiments have shown that 
the order parameter of YBCO does indeed have a signif- 
icant s-wave component [sof . 

In summary, we have presented a group theoretic anal- 
ysis of the several organic superconductors. This analysis 
has lead us to propose order parameters summarised in 
table Uni In detail we have argued that orthorhombic 
^LOS (e.g., K-Br and k-C1) have 'c?a;2_j,2' order parame- 
ters but that monoclinic ^LOS (e.g., k-NCS) have nodes 
along the c-axis and possibly another node in the con- 
ducting plane that is not required by symmetry. We pro- 
posed that highly frustrated ^LOS, such as K-CN3 and 
/?'-dmit, have ^dx2_y2 -\- idxy order parameters but the 
transition is split by the symmetry of the crystal so that 
there is first a transition from the normal state to a 'd- 
wave' superconductor and then a second phase transition 
from the 'd-wave' superconductor to a 'c? -f id' supercon- 
ductor. We argued that ^LOS, e.g., 9-lz, have 'd^y -I- s' 
order parameters. 
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